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We study the steady state structure and dynamics of an interface in a pure Ising system on a square 
lattice placed in an inhomogeneous external field with a profile designed to stabilize a flat interface, 
and translated with velocity We. For small We, the interface is stuck to the profile, is macroscopically 
smooth, and is rippled with a periodicity in general incommensurate with the lattice parameter. 
For arbitrary orientations of the profile, the local slope of the interface locks in to one of infinitely 
many rational values (devil's staircase) which most closely approximates the profile. These "lock-in" 
structures and ripples dissappear as We increases. For still larger We the profile detaches from the 
interface. 

PACS numbers: 05.10Gg,64.60.Ht,68.35.Rh 



The ability to grow flat solid surfaces [l| is often of ma- 
jor technological concern, for example, in the fabrication 
of magnetic materials for recording devices where sur- 
face roughness[2] causes a sharp deterioration of mag- 
netic properties. Most growing surfaces or interfaces on 
the other hand coarsen[3j with a width, cr, which diverges 
with system size as L" and time as , where a and j3 
(a, /? > 0) are the roughness and dynamical exponents 
respectively. Is it possible to drive such an interface with 
a pre-determined velocity vj and, at the same time, keep 
it flat (i.e. a = /3 = 0) ? In this Letter, we explore 
this possiblity by studying growing interfaces in a non- 
uniform field with an appropriately shaped profile mov- 
ing without change of shape at an externally controllable 
velocity Ve- We find that for less than a limiting value 
Woo, it is possible to produce a macroscopically flat, in- 
terface growing with average velocity Vf = Vg. Micro- 
scopically, however, the interface shows an infinity of dy- 
namical ripple structures similar to self-assembled, com- 
mensurate -incommensurate(C-I) [1] domains produced 
by atomic mismatch [B']. The ripples vanish with increas- 
ing Ve through a fluctuation induced C-I transition. 

Real solid-solid interfaces being complex [6|, an under- 
standing of the dynamics of such interfaces in a gen- 
eral time -dependent, inhomogeneous field may only be 
achieved by beginning with a relatively simple, but non- 
trivial, system viz. an interface in the two -dimensional 
(2-d) Ising model. Our results, therefore, concern mainly 
this model system, though towards the end we discuss 
possible modifications, if any, for solid interfaces. 

The dynamics of an Ising interface in a (square) lat- 
tice driven by uniform external fields is a rather well 
studiedjl] subject. The velocity, Voo of the interface 
depends on the applied field, h and the orientation 9 
measured with respect to the underlying lattice. The 
interface is rough and coarsens with KPZ(3 exponents 
a = 1/2 and (3 = 1/3. 

Consider, an interface Y{x, t) between phases with 



magnetizatiori, (/)(a;, y, t) > and 4>{x,y,t) < 0, in a 2-d 
square lattice [7| obeying single-spin flip Glauber dynam- 
ics in the limit h/J,T/J 0. Here J is the Ising 
exchange coupling and T the temperature. An exter- 
nal non-uniform field with a profile h{y,t) = hmaxfiy,t) 
where f{y,t) = tanh((?/ — Vf,t)/x) and x is the width of 
the profile (see Fig[T](a)) is applied such that h = hmax in 
the > and —hmax in the 4> < regions separated by 
a sharp edge. The driving force depends on the relative 
local position of Y(x, t) and the edge. In the low temper- 
ature limit the interface moves solely by random corner 
flips 3] (Fig. [Ijb)), the fluctuations necessary for nucle- 
ating islands of the minority phase in any region being 
absent. We study the behaviour of vj and the structure 
of the interface as a function of Ve and orientation. 

Naively, one expects fluctuations of the interfacial co- 
ordinate Y{x, t) to be completely suppressed in the pres- 
ence of h{y,t). Indeed, as we show below, a mean field 
theory obtains the exact behaviour oivf as a function of 
Ve (Fig. [11(c)). Using model A dynamics^ for a coarse- 
grained Hamiltonian (which, for the moment, ignores the 
lattice) of an Ising system in a external non-uniform field 
together with the assumption that the magnetization (p 
is uniform everywhere except near the interface, one can 
derive[ll| an equation of motion for the interface. 



^ - A — 

dt dx"^ 



'A2(|^)'/(r,i)-A3/(r,i)+c(x,t) (1) 



where Ai,A2 and A3 are parameters and C is a Gaus- 
sian white noise with zero mean. Note that Eq. ^ lacks 
Galilean invariance[lo| Y' Y+ Azet, t' t. 

A mean field calculation amounts to taking Y = Y{t) 
i.e. neglecting spatial fluctuations of the interface and 
noise. For large times (i 00), Y — > Vft, where 
Vf is obtained by solving the self-consistency equation; 
Vf = limt_oo -A3tanh{(w/-Ve)t/x} = -A3 sign(w/-We). 
For small Ve the only solution (Fig.[ljc), inset) is Vf — Ve 
and for Ve > Voo, where Uqo — A3 we get Vf ~ X3 ~ Vrx,. 
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FIG. 1: (a) An Ising interface Y{x,t) (bold curved line) be- 
tween regions of positive (marked +) and negative (marked 
— ) magnetization in an external, inhomogeneous field with 
a profile which is as shown(dashed line). The positions of 
the edge of the field profile and that of the interface are la- 
belled Se and Sf respectively, (b) A portion of the interface 
in a square lattice showing a corner, (c) The interface ve- 
locity w/ as a function of the velocity of the dragging edge 
Ve for Ns = 100(a), 1000(0), lOOOO(-h) and p = 0.5. All the 
data (0,0,-1-) collapse on the mean field solution (dashed 
line). Inset shows the graphical solution (circled) of the self- 
consistency equation for vf, dashed line represents Vf = Vf. 



We thus have a sharp transition (Fig.[TJc)) from a region 
where the interface is stuck to the edge to one where it 
moves with a constant velocity. How is this result altered 
by including spatial fluctuations of Y 7 This question is 
best answered by mapping the interface problem to a 1-d 
cellular automaton [3l Il2|. 

The interface coordinate Y{x,t) in a square lattice is 
represented [l^ [isj by the set of integers {yi}, I < i < 
Np denoting positions of Np hard-core {ui+i > Vi + I) 
particles in a 1-d lattice of Ns sites. The particle density 
p = Np/Ns determines the mean slope of the interface 
tan Of = l/p and the motion of the interface by corner 
flips corresponds to the hopping of particles with right 
and left jump probabilities p and q (p + q = 1). Trial 
moves are attempted sequentially on randomly chosen 
particles [l3| and Np attempted hops constitute a single 
time step. In our case, p and q are position dependent 
such that, Ai{ve t) ^ p — q — A sign(yi — i/ p — Vet) with 
A = 1. The bias Ai{vet) is appropriate for a step func- 
tion (x = 0) proflle with the slope of the profile edge 



equal to the average slope of the interface. We track 
the instantaneous particle velocity Vf{t) defined as the 
number of particles moving right per unit time, the av- 
erage position < y{t) >= Np^J2i=iN yi{^)^ width 

a2(i) = A^-iEz=i,jv, < {y^{t)- < y^ih >? > and the 
local slope of the interface given by the local density of 
particles. Angular brackets denotes an average over the 
realizations of the random noise. The usu al p article hole 
symmetry for an exclusion process d, [H, 15 1 is violated 
since exchanging particles and holes alters the relative 
position of the interface compared to the edge. 

We perform numerical simulations [l3| of the above 
model for Ng upto lO"' to obtain Vf for the steady state 
interface as a function of Ve as shown in Fig. [Ijc). A 
sharp dynamical transition from an initially stuck in- 
terface with Vf — to & free, detached interface with 
Vf = Vqo = A{1 — p) is clearly evident as predicted by 
mean field theory. The detached interface coarsens with 
KPZ exponents [11|. Note that, even though the mean 
field solution for Vf{ve) neglects the fluctuations present 
in our simulation, it is exact. The detailed nature of 
the stuck phase (w/ — Ve and a bounded) is, on the other 
hand, considerably more complicated than the mean field 
assumption Y{x^t) — Y{t). 

The ground state of the interface in the presence of 
a stationary (ug = 0) field profile is obtained by mini- 
mizing E — ^/NpJ2i{yi ~ i/ P ~ cY with respect to the 
set {yi\ and the constant, c, subject to the constraint 
that yi are integers. The form of E leads to an ad- 
ditional non-local, repulsive, interaction between parti- 
cles. The minimized energy may be calculated exactly, 

E{p — m/n) — i(i — — ^) + t'^ ir-^ for m even 

and ]i^(l — for ni odd, where the density p = m/n, 
is a rational fraction. The energy satisfies the bounds 
E{l/n) = < E{m/n) < limm^oo E (m/n) = 1/12 
where the upper bound is for an irrational density. For an 
arbitrary < p < 1 the system ({yi}) prefers to distort, 
conforming within local regions, to the nearest low- lying 
rational slope l/p interspersed with an ordered array of 
"discommensurations" of density pd = \p — p\ which ap- 
pear as long wavelength ripples (see Fig.[3]inset (c)). The 
p as a function of p shows a "devil's staircase" structure 
(complete for 0) with a multifractal [l3| measure. 

We observe this in our simulation by analyzing the in- 
stantaneous distribution of the local density of particles 
to obtain weights for various simple rational fractions 
upto generation g = 9 in the Farey tree of rationals 14 1 . 
A time average of the density corresponding to the frac- 
tion with largest weight at any t, then give us the most 
probable density p — distinct from the average p which 
is constrained to be the inverse slope of the profile. For 
small Ve the interface is more or less locked in to a single 
p, shown as white regions in the phase diagram (Fig. [2]) 
in the v^ — p plane. 

For low velocities and density where correlation effects 
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FIG. 2: The dynamical phase diagram in and p plane. The 
numbers on the p axis mark the fractions p, which determines 
the orientation of the lock-in phase. The three regions white, 
black and grey correspond to the rippled, the disordered and 
the detached phases respectively. 



due to the hard core constraint are negligible, the dynam- 
ics of the interface may be obtained exactly[llj. Under 
these circumstances the Np particle probability distribu- 
tion for the yi's, P{yi,y2, ■ ■ ■ ,yNp) factorizes into single 
particle terms P{yi). Knowing the time development of 
P{yi) and the ground state structure the motion of the 
interface at subsequent times may be trivially computed 
as a sum of single particle motions. A single particle 
(with say index i) moves with the bias Ai(wei) which, 
in general, may change sign at y < i/p + Vet < y + 1. 



Solving the appropriate set of master equations we 
obtain, for << 1 the rather obvious steady state solu- 
tionP(yi) — l/2((5j^._y-|-(5j^.^2^+i) and the particle oscillates 
between ?/ and y -1-1. Subsequently, when i/p-f-Ugi > 
the particle jumps to the next position and P{yi) relaxes 
exponentially with a time constant t = 1 to it's new 
value with y ^ y -\- 1. In general, for rational p = m/n, 
the motion of the interface is composed of the indepen- 
dent motions of m particles each separated by a time lag 
of Ti = \/mVe- The result of the analytic calculation 
for small Ve and p has been compared to those from sim- 
ulations in Fig. [3] for p = 1/5 and 2/5. For a general 
irrational p < 1/2, m ^ oo consequently, ^ 0. The 
yiS are distributed uniformly around the mean imply- 
ing (T^ = 1/3 independent of system size and time. For 
p > 1/2 the width cr^ = (1 — p)/3p since the number of 
mobile particles decreases by a factor of (1 — p)/p. 

The forward motion of an irrational interface is accom- 
panied by the motion of discommensurations along the 
interface with velocity Ve which constitutes a kinematic 
wave0, [3] parallel to the interface. As the velocity Ve 
is increased the system finds it increasingly difficult to 
maintain its ground state structure and for t > tl the in- 
stantaneous value of p begins to make excursions to other 



FIG. 3: Variation of < y{t) > with t for = 0.025 and 
p = 1.0. Lines denote analytic results while points denote 
Monte Carlo data for p = 1/5 (uppermost curve), 2/5 and 
an incommensurate p near 1/3. Inset (a)-(c) shows the cor- 
responding ground state interfaces (j/i — i/p). The arrows in 
(c) mark the positions of two discommensurations. 



nearby low-lying fractions and eventually becomes free. 
Steps corresponding to p = m/n dissappear (i.e. p ^ p) 
sequentially in order of decreasing m and the interface 
loses the ripples. The interface is disordered though a 
and (3 continue to be zero (black region in Fig. ^ . The 
locus of the discontinuities (within an accuracy of 1/iVp) 
in the p(p) curve for various velocities Ve gives the limit 
of stability of the lock-in rippled phases. 

While the stability of mismatch domains Q is decided, 
mainly, by competition between mechanical, long-ranged 
(elastic) and short-ranged (atomistic) interactions dy- 
namical ripples vanish with increasing Ve through in- 
creased fluctuations. We argue that it is sufficient to 
project the entire configuration space yi of the stuck in- 
terface onto the single variable p. In Fig. 3] (inset) we 
have plotted the energy E{p) for the ground-state con- 
figuration with density p. As is obvious from the figure, 
E{p) has a structure similar to the free energy surface 
of a 1-d "trap" model[l3| often used to describe glassy 
dynamics. The distinction, of course, is the fact that 
the energies of the traps in this case are highly corre- 
lated. We may then describe the dynamics of the stuck 
interface as the Langevin dynamics [9| of a single particle 
with coordinate p' diffusing on a energy surface given by, 
F{p') = E{p') + k{p' — p)^, kicked by a Gaussian white 
noise of strength T. The second term, containing the 
modulus K, ensures that p' — > p for i ^ oo. At interme- 
diate times, however, the system may get trapped indef- 
initely in some nearby low-lying minimum with p' — p 
if T (oc Wg by symmetry) is not large enough. As T in- 
creases, the time spent in jumping between minima may 
exceed the residence time in the minimum resulting in a 
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FIG. 4: A surface plot of p vs p and VTk showing steps 
for the fractions 1/4,2/7 and 1/3. Note that the step corre- 
sponding to 2/7 vanishes at \/Tk > 0.5. Inset (top) shows a 
plot of E{p) for a, Ns = 420 site system. 



noise induced C-I transition (Fig. d]) from p. 

In this Letter, we have studied the static and dynamic 
properties of an Ising interface in 2-d subject to a non- 
uniform, time-dependent external magnetic field. The 
system has a rich dynamical phase diagram with in- 
finitely many steady states. The nature of these steady 
states and their detailed dynamics depend on the orien- 
tation of the interface and the velocity of the external 
field profile. How are our results expected to change 
for real driven solid interfaces? Firstly, real field pro- 
files would have a finite width x > 0. However, as long 
as X is comparable to atomic dimensions we find no ap- 
preciable change in the results. Secondly, elastic interac- 
tion between "particles" or steps in the interface [6] may 
smoothen the devil's staircase though we expect that for 
typical solids this will be minimal. Thirdly, the struc- 
ture of the underlying lattice may influence the stability 
of particular orientations. Finally, exciting new physics 
may come into play as new modes e.g. point and line 
defects [l7[, as well as phonon degrees of freedom (leading 
to acoustic emmissions 6^ ) are accessed as v,, increases. 
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